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PREFACE 


The  work  described  in  this  report  was  authorized  under  Project  No.  0602622A552. 
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official  endorsement  of  any  commercial  products.  This  report  may  not  be  cited  for  purposes  of 
advertisement. 


This  report  has  been  approved  for  public  release.  Registered  users  should  request 
additional  copies  from  the  Defense  Technical  Information  Center;  unregistered  users  should  direct  such 
requests  to  the  National  Technical  Information  Center. 
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EMISSIVE  VERSUS  ATTENUATING  SMOKES 


1.  INTRODUCTION 

Theory  is  developed  and  one  dimensional  radiative  transfer  calcailations  are  carried  out 
predicting  contrast  and  signal  to  noise  transmittances  of  two  types  of  smoke  with  the  purpose  of 
determining  what  reductions  in  optical  depth  would  be  possible  when  smoke  thermal  emission  is 
increased.  First,  analysis  was  convicted  on  scattering,  absorbing  and  emitting  smoke  clouds  at  elevated 
temperatures  where  entrained  air  terrq)erature  is  equal  to  aerosol  particle  temperature.  Second,  analysis 
was  completed  on  smoke  clouds  containing  both  scattering  nonabsorbing  smoke  particles  at  ambient 
temperature  and  high  temperature  emitting  flare  particles.  In  the  first  case  it  was  found  that  heating  the 
smoke  particles  will  not  significantly  reduce  contrast  transmittance  or  signal  to  noise  ratio  transmittance 
because  smoke  particles  dump  heat  into  the  air  within  a  few  microseconds  and  even  the  most  exothermic 
reactions  are  not  energetic  enough  to  heat  the  entrained  air  mass  more  than  a  few  degrees  Kelvin  at  smoke 
battlefield  concentrations  of  a  few  tenths  of  a  gram  of  aerosol  per  cubic  meter  of  air.  If  the  heat  comes 
from  another  source  such  as  smoke  particles  injected  into  a  hot  stream  of  engine  exhaust  gas  or 
pyrotechnically  produced  gas,  the  cloud  temperature  would  have  to  be  elevated  to  a  level  where  it  would 
become  buoyant  and  rise  rapidly  before  significant  reductions  in  optical  depth  could  be  realized.  In  the 
second  case,  if  flare  particles  are  much  larger  than  the  smoke  particles  they  will  dissipate  heat  primarily 
by  radiation  rather  than  conduction  into  air  and  much  higher  ten:5)eratures  and  thermal  emission  can  be 
achieved.  On  the  other  hand  when  flare  particles  are  made  larger  their  emission  cross  section  per  mass 
decreases  and  the  gravitational  settling  velocity  increases.  Therefore  flare  particle  size,  shape  and 
reaction  rate  or  tenqjerature  must  be  optimized  to  balance  maximization  of  flare  particle  emission  and 
duration  with  minimization  of  attenuating  particle  cloud  optical  depth  maintaining  contrast  and  signal  to 
noise  ratio  transmittance  at  an  acceptable  level.  The  optical  depth  requirement  can  be  interpreted  as  a 
requirement  on  either  extinction  coefficient  or  concentration  pathlength  product  since  the  optical  depth  is 
the  product  of  both.  Calculations  based  on  one  dimensional  radiative  transfer  are  most  accurate  for 
extensive  smoke  blankets  and  walls  where  photon  flux  passing  through  cloud  edges  can  be  neglected. 


2.  CONTRAST 

The  ability  to  see  an  object  depends  on  the  contrast  between  the  object  and  backgroxmd.* 
When  smoke  is  placed  between  the  observer  and  object  the  contrast  is  affected  in  two  ways;  first  the 
smoke  attenuates  the  radiance  of  the  object  and  backgroxmd,  second  the  smoke  superinposes  scattered 
and  emitted  radiance  onto  the  scene.  Smoke  attenuation  of  background  and  object  radiance  is  predicted 
by  the  smoke's  unscattered  transmittance  while  radiance  superinposed  onto  the  scene  is  predicted  by 
summing  contributions  due  to  reflected,  scattered  transmitted  and  emitted  radiance.  In  order  to  quantify 
the  effects  of  attenuation  vs  scatter  plus  emission  on  contrast  we  first  recall  a  standard  definition  of 
contrast  of  an  object  against  a  background^ 


C=(L„-L,)/(L„  +  L,) 


where  Lo  is  the  object  radiance  and  Lb  his  background  radiant  and  then  write  the  contrast  of  that  object 
against  that  background  viewed  through  smoke 

c'=(l;-l;)/(l;+l;) 

where  the  object  and  background  radiances  viewed  through  smoke  are 


7 


L„'  =  T„L„  +  T,L^  +  RL_+Eo-T"  and  lJ  =  T„Lb +T3L^ +RL_ +  EctT'‘ 


Here  Tu  is  the  unscattered  (Beer's  Law)  transmittance,  T,  is  the  scattered  light 
transmittance,  R  is  the  reflectance,  L_  is  ambient  sky  plus  solar  phis  background  terrain  radiance  directed 
away  from  the  observer,  L+  is  ambient  sky  plus  solar  plus  background  terrain  radiance  directed  toward  the 
observer,  E  is  the  smoke  emitted  radiance,  cr  is  the  Stefan  Bolt2man  constant  and  T  is  the  cloud 
temperature.  The  unscattered  transmittance  through  the  smoke  depends  only  on  the  extinction  coefficient, 
extinction  cross  section  per  mass  of  smoke,  and  the  concentration  path  length  product.  The  smoke 
reflectance  and  scattered  light  transmittance  depend  on  these  parameters  as  well  as  single  scatter  albedo, 
the  ratio  of  scatter  cross  section  to  extinction  cross  section  and  phase  function  (the  differential  scattering 
pattern).  Smoke  emission  depends  on  all  of  the  previously  mentioned  parameters  as  well  as  the  cloud 
temperature  in  the  case  of  thermal  emission  and  other  properties  in  the  case  of  fluorescence. 


3.  CONTRAST  TRANSMTTTANCE  AND  ONE  DIMENSIONAL  RADIATIVE 

TRANSFER  MODEL 


To  quickly  analyze  the  effects  of  all  these  smoke  properties  upon  contrast  it  is  convenient 
to  work  with  a  one  dimensional  radiative  transfer  model.  Here  the  smoke  is  considered  as  a  smoke  wall 
or  smoke  blanket  and  the  analysis  requires  considerably  less  computation  time  due  to  the  single  closed 
form  nature  of  these  solutions  rather  than  the  numerical  solutions  involved  with  three-dimensional 
radiative  transfer,  bn  this  way  we  dispense  with  considerations  of  cloud  shape  and  transport  of  radiation 
perpendicular  to  the  line  of  sight,  non  uniform  backgrounds  and  non  uniform  ambient  illumination. 
Further  sin^rlifications  that  speed  computation  involve  replacing  the  phase  function  with  a  single  number, 
the  asymmetry  factor  and  reducing  the  flow  of  radiation  to  diffuse  and  collimated  conqjonents  moving  in 
both  directions  by  using  the  four  flux  model  of  one-dimensional  radiative  transfer  where  the  flux  is 
considered  as  a  combination  of  incident  coUimated  flux  passing  through  the  smoke  layer  in  both 
directions  along  with  scattered  flux  going  in  both  directions.  Diffuse  incident  flux  can  be  represented  by 
averaging  collimated  flux  over  angle. 

Two  simplifications  result  when  we  work  with  the  contrast  transmittance  Ct,  ratio  of 
contrast  observed  through  smoke  over  contrast  observed  without  smoke  present, 

Lq  +  lQ  1 

L„-4  1,  ,  2T,U+2RL.+2EaT^ 

T„(L„+L,) 

rather  than  just  contrast.  First  the  range  of  Cj  always  lies  between  0  and  1 .  Second  the  set  of  four 
radiance  ratios  Lb/Lo,  L+/Lo,  L  VLo  and  E/Lq  reduce  to  three;  L+/(Lo+Lb),  L  7(Lo+Lb)  and  E/(Lo+Lb).  The 
four  flux  theory  one  dimensional  radiative  transfer  theory  connects  the  smoke  blanket  or  smoke  wall 
imscattered  light  transmittance  Tu,  reflectance  R,  scattered  light  transmittance  Ts  and  emittance  E  with 
smoke  particle  properties,  aerosol  concentration  and  cloud  pathlength^ 


Cx  = 


T..(L„-LJ 


T.(L +L.)  +  2XL  +2RL  +2EcrX 


R  = 


(1-p  +  z) 
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where 


T  = 


2  z 


(l-p  +  zje%  -  ^l-p-z^e-% 
E  =  (l-R-T,-T„) 


-Y 

.g-A 


p=|0+g) 

q^O-g) 

z  =  ((l-a)(l-ag))^. 

The  smoke  optical  dq)th  is  designated  as  b,  single  scatter  albedo  is  designated  as  a,  asymmetry  factor  is  g, 
and  jU  is  the  cosine  of  the  angle  of  incidence  with  respect  to  the  cloud  wall  or  blanket  surface. 


4.  SIGNAL  TO  NOISE  RATIO  AND  SIGNAL  TO  NOISE  TRANSMITTANCE 

For  photon  noise  limited  detection  a  signal  to  noise  ratio  SN  is  defined  in  terms  of  the 
number  of  photons  collected  n  which  is  proportional  to  the  collected  radiance  L 


n  =  /?L 


where  is  the  proportionality  constant  incorporating  the  collection  area  and  time.**’^  Setting  the  noise 
equal  to  the  square  root  of  the  number  of  photons  collected  in  each  measurement  and  the  total  noise  of 
object  and  background  measurements  equal  to  the  square  root  of  the  sum  of  the  squares  we  have 


(t„L„  +  T,L,  +  RL.  +Eorr ) - (t„L,  +T,L,  +  RL_  +  ErrT^) 

SN  =  -X - L 

^y5[(T„L„  +  T.L,  +RL_  +Ec7r  )+(T„L,  +  T.L,  +  RL.  +  E(Tr  )] 

=  (L„  -L,)/  (T„(L„+L,)+2T,L,  +  2RL_  +2E^tT^)^ 


We  define  signal  to  noise  transmittance  SNt  as  the  ratio  of  signal  to  noise  ratio  looking  through  smoke 
divided  by  the  signal  to  noise  ratio  without  smoke  present. 


.  V^T,(L,-L,) _ (L„+L,) 

"  (t„(L„  +L,)  +  2T,L,  +2RL_  +2E<Tr  ) 
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1 


SN,  =  ■ 


rK 


1  + 


2TL 


2RL 


+  - 


2EaT 


4  Vi  f 

'^TT  I  noT  .  ^y2 


t„(l„+4)  t„(l„  +  l,)  t„(l„  +  l,) 


1+ 


2TL  +2RL_  +  2Ecrr 


T„(k  +  L,) 


Contrast  transmittance  is  simply  related  to  signal  to  noise  transmittance 

SN,=(T„C,)>^ 


Three  simplifications  result  when  we  work  with  signal  to  noise  transmittance  instead  of 
the  signal  to  noise  ratio.  First,  the  range  of  SNt  lies  between  0  and  1 .  Second,  the  set  of  four  radiance 
ratios  Li/Lo,  L  JLq,  L_/Lo  and  E/Lo  reduce  to  three;  L+/(Lo+Lb),  L_/(Lo+Lb)  and  E/(Lo+Lb).  Finally,  the 
proportionality  constant  y9can  be  eliminated. 

It  is  clear  from  these  expressions  that  contrast  transmittance  and  signal  to  noise 
transmittance  are  reduced  by  increasing  cloud  reflectance  and  scattered  light  transmittance  but  they  are 
also  reduced  by  increasing  cloud  emittance  and  temperature.  Two  parameters  are  defined  to  quantify 
decreases  in  contrast  and  signal  to  noise  transmittance  due  to  heating  the  cloud  relative  to  those  decreases 
produced  by  an  identical  unheated  cloud.  These  parameters  are  the  contrast  transmittance  ratio  CTR  and 
the  signal  to  noise  transmittance  ratio  SNTR;  the  heated  cloud  value  divided  by  the  unheated  cloud  value. 


5.  HEATING  SMOKE  PARTICLES  BY  REACTIONS  AND  EXTERNAL  MEANS 

Heating  smoke  particles  by  an  exothermic  reaction  or  using  external  radiation  from  an 
energetic  beam  such  as  a  high  energy  laser  will  quickly  lead  to  heat  loss  to  the  air  mass  within  the  smoke 
cloud  which  is  generally  about  four  orders  of  magnitude  greater  than  the  aerosol  mass  within  the  cloud. 
Exothermic  smoke  reactions  producing  AH  joules/gram  of  smoke  particles  will  quickly  elevate  smoke 
cloud  temperature  by  AT  by 


AT  =  AHC/(pc) 

where  C  is  smoke  concentration  in  grams/m^  p  is  air  density  (~ 1,300  g/m^)  and  c  is  the  heat  capacity  of 
air  (~1  j/g°K).  In  order  to  raise  smoke  cloud  temperatiare  by  more  than  AT  a  reaction  must  produce  AH 
greater  than 


AH  =  /9cAT/C 

For  example,  a  smoke  particle  reaction  would  have  to  produce  greater  than  20,000  J/g  of 
heat  to  raise  cloud  temperature  more  than  3°K  at  smoke  concentrations  under  0.2  g/m^.  Similarly  we  can 
heat  the  aerosol  cloud  using  an  external  source  of  radiation  such  as  a  high  energy  laser  or  spotlight.  The 
temperature  rise  depends  on  the  absorbed  radiance  in  accordance  with  the  energy  balance  equation 

TuE  %C  Adl  =  cp  ATAdl 

where  Tu  is  the  unscattered  light  transmittance,  E  is  the  incident  high  energy  density  in  joules/m^  %is  the 
aerosol  absorption  coefficient  (absorption  cross  section  per  aerosol  mass),  A  is  the  cloud  area  transverse 
to  radiance  propagation  and  dl  is  the  incremental  pathlength  in  the  cloud  along  the  beam.  Thus  to  raise 
smoke  temperature  by  more  than  AT  the  high  energy  density  must  exceed 
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E  >  c  /?  AT  /  (  Tu  C  ) 


For  example,  incident  energy  densities  greater  than  19,500  J/m^,  achievable  with  a  one 
meter  diameter  19.5  kilowatt  beam  incident  on  a  cloud  in  a  1  meter  per  second  cross  wind,  would  be 
required  to  raise  the  temperature  of  the  cloud  edge  facing  the  incident  power  by  3"K  for  a  cloud  with  an 
absorption  coefficient  of  1  mVg  at  a  concentration  of  0.2  g/m^.  Finally  a  third  approach  would  involve 
heating  the  aerosol  and  air  by  mixing  with  a  hot  exhaust  gas  stream  or  pyrotechnically  generated  hot  gas 
flow.  Now  the  external  power  P  that  is  required  to  elevate  the  temperature  by  AT  is  given  by 


P  =  cihAT 


where  rii  is  the  mass  flow  rate  of  hot  exhaust  or  pyrotechnically  generated  gas  plus  aerosol  plus  entrained 
ambient  air  that  is  mixed  producing  the  smoke  cloud.  For  example,  6  kilowatts  of  power  must  be 
provided  to  elevate  a  mass  flow  rate  of  2  kilograms  per  second  by  3  °K. 


6.  HEATED  SMOKE  CLOUD  BUOYANCY 

Heated  clouds  will  of  course  rise  due  to  their  buoyancy.  The  velocity  of  rise  can  be 
estimated  by  balancing  the  buoyancy  force  acting  on  the  cloud 

FB  =  A^(7T/6)D^g 


with  the  drag  force  acting  on  the  cloud 


Fd  =  Cd(7t/4)D^/?v^/2 

where  is  the  reduction  in  cloud  density  due  to  heating,  D  is  the  cloud  diameter,  g  is  acceleration  due  to 
gravity  (~10m/s^),  Cd  is  the  air  drag  coefficient  acting  on  the  entire  cloud  and  v  is  the  velocity  that  the 
cloud  is  rising.  The  ideal  gas  law  predicts  the  reduction  in  cloud  density  due  to  a  temperature  elevation 
AT  above  ambient  temperature  T 


V=/?/(1+T/AT) 

The  drag  coefficient  for  Reynolds  numbers  greater  than  200,000  is  0.1 .  A  correction 
factor  that  reduces  this  drag  coefficient  due  to  internal  air  flow  within  the  cloud  can  be  computed  in  the 
laminar  flow  regime  and  gives  the  value  5/6  for  a  gas  sphere  moving  through  a  gas  of  equal  viscosity.* 
Although  we  are  outside  the  laminar  flow  regime  we  can  expect  the  cloud  to  break  up  anyway  once 
velocities  become  great  enough  to  produce  significant  turbulence.  The  Reynolds  number  Re  is  given  by 

Re  =  V  D  /  // 

where  n  is  the  kinematic  viscosity  of  air  (0.152  cmVs  at  293°K  ).  Setting  the  buoyancy  force  equal  to  the 
air  drag  force  and  solving  for  the  velocity  of  cloud  rise  we  have 

v  =  ((Ap(7r/6)D^g)/(CD(;T/4)D'p/2)y''  =(4Ap  gD/(3pCD)y'' 

Substituting  for  Lp  we  get 
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For  example,  a  10  meter  diameter  cloud  with  a  temperature  elevation  of  3"K  above  an  ambient 
temperature  of  300”K  would  rise  at  a  velocity  of  3.6  m/s.  This  corresponds  to  a  Reynolds  number 

Re  =  (3.6  m/s)  (10  m)  /  (0.0000152  m^/s)  ~  2,370,000 

which  is  in  the  range  where  the  drag  coefficient  is  approximately  equal  to  0.1  but  we  can  expect  turbulent 
flow  to  quickly  break  apart  the  cloud  in  this  flow  regime. 


7.  TIME  SCALE  FOR  HEAT  TRANSFER  FROM  AEROSOL  PARTICLES 


The  rate  of  heat  conduction  away  from  the  smoke  particles  into  the  air  surrounding  the 
particles  can  be  estimated  using  the  instantaneous  or  continuous  point  heat  source  solutions  for  transient 
heat  conduction  into  an  infinite  medium.  The  instantaneous  point  source  solution  is’ 


AT(r,t)=- 


8(mrt)^/3c 


O 

V  -  Aki 


where  A- is  the  thermal  difihisivity  of  air,  0.187  cm’/s,  Q  is  the  quantity  of  heat  instantaneously  released  at 
time  t  =  0  and  r  is  the  distance  from  the  point  source.  The  temperature  has  its  maximum  value  at  t=r^/(6A-) 
which  is  also  the  time  at  which  the  heat  has  reached  a  mean  square  distance  r^  from  the  instantaneous 
source.  The  continuous  point  source  solution  is’ 


AT(r,t)  =  — ^ — etfc  S — 

^TTKTpc  v4/rt 

where  Q  is  the  rate  of  release  by  the  continuous  point  source  and  etfc  is  the  complimentary  error  function 
which  can  be  approximated  by  the  expression* 


I  ^ 

etfc(x)=  1- Vl-e 

After  a  time  t  =  r’/(6A')  the  heat  initially  released  by  the  continuous  source  has  moved  a 
mean  square  distance  r’  away  from  the  source  and  etfc{(3/2f'^)  =  0.077,  which  is  still  a  small  fraction  of 
the  steady  state  solution 


AT(r,  t  =oo)  =  - - 

Ammpc 

It  is  important  to  recognize  that  the  instantaneous,  continuous  and  steady  state  point 
source  solutions  apply  to  ensembles  of  particles  up  to  times  where  there  is  insignificant  heating  of  an  air 
parcel  by  more  distant  particles  relative  to  heating  by  the  nearest  particle.  In  order  to  determine  how 
quickly  an  aerosol  particle  dumps  heat  into  the  air  we  compute  the  time  required  for  a  particle  to  heat  air 
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within  the  cloud  located  a  maximum  distance  from  any  smoke  particle.  Therefore  we  must  estimate 
interparticle  spacing  within  a  cloud.  For  typical  smoke  clouds  we  assume  a  particle  volume  of  one  cubic 
micrometer  and  a  typical  concentration  of  0.2g/m^  If  we  assume  a  typical  particle  density  of  2  then  each 
particle  will  weigh  2X10'’^g  and  there  will  be  lO”  particles  in  a  cubic  meter  of  aerosol.  If  the 
concentration  is  homogeneous  the  spacing  between  particles  will  be  about  10’**^^  meter  or  214  |nn  and 
maximum  distances  from  any  given  particle  will  be  about  half  that  value.  The  maximum  time  required 
for  heat  to  travel  to  any  parcel  of  air  within  the  cloud  will  correspond  to  a  mean  square  distance  equal  to 
the  maximum  distance  away  from  any  particle  in  the  cloud,  107  ^im, 

t=(0.0107cm)V((6)(0.187cmVs))=  102  psec. 

Thermal  conduction  equilibrium  is  achieved  within  several  multiples  of  this  time,  probably  less  than  one 
millisecond. 


8.  FLARE  PARTICLE  HEAT  TRANSFER 

If  we  want  to  partition  more  heat  lose  into  radiation  that  can  reduce  contrast  and  signal  to 
noise  transmittance  rather  than  conduction  into  the  air,  which  has  no  effect  on  contrast  or  signal  to  noise 
ratio  in  atmospheric  windows,  we  can  increase  particle  size  and  temperature.  Increasing  particle  size  will 
decrease  surface  area  available  for  both  conductive  and  radiative  losses  and  increase  the  duration  of 
particle  reactions  and  increase  settling  velocity.  To  see  how  particle  size  influences  particle  temperature 
we  first  look  at  the  steady  state  solution  for  conduction  of  heat,  for  example  generated  by  surface  and 
volume  reactions  in  a  spherical  particle  of  diameter  f ,  into  the  surrounding  air.  We  assume  that  the 
spacing  between  these  large  “flare”  particles  is  great  enough  that  heat  conduction  by  any  given  flare 
particle  is  imaffected  by  heat  conduction  by  surrounding  particles.  The  temperature  elevation  of  the 
surface  of  a  spherical  particle  producing  AH^  watts  of  power  per  unit  particle  volume  will  be 
proportional  to  the  particle  size  squared’ 


AT,=AH,#’/3K 

while  for  surface  heating  producing  AH^  watts  of  power  per  imit  particle  surface  area  the  temperature 
elevation  is  linearly  proportional  to  particle  size’ 

AT3=AH,^/K 

where  K=x/7  c  is  the  conductivity  of  air.  Volume  heating  could  be  the  result  of  either  exothermic 
chemical  reactions  occurring  throughout  the  particle  volume  or  absorption  of  incident  high  power  density 
radiation  throughout  the  particle  volume,  which  is  often  the  case  when  particle  size  is  small  compared  to 
the  wavelength  of  the  high  power  radiation.  Surface  heating  could  be  the  result  of  either  exothermic 
chemical  reactions  occurring  on  the  particle  surface  or  absorption  of  incident  high  power  density  radiation 
over  the  particle  surface,  which  is  the  case  when  particle  size  is  large  compared  to  the  wavelength  of  the 
high  power  radiation.  If  we  now  assume  all  losses  are  radiative  rather  than  conductive,  energy  balance 
predicts  the  equilibrium  temperature  for  volume  heating  will  be  proportional  to 

T.  =(aH,#/(3  <rs)f 

and  the  equilibrium  temperature  for  surface  heating  will  be  independent  of  particle  size 
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T,  =  (AH7Wr 

where  e  is  the  emissivity  or  absorptivity  of  the  particle.  These  two  expressions  for  emitted  radiation 
show  clearly  the  particle  size  dependence  of  equilibrium  temperatures  for  particles  that  radiate  from  their 
surface  rather  than  their  volume.  This  is  true  for  particles  larger  than  the  wavelength  of  emitted  radiation 
because  the  emissivity  does  not  depend  on  particle  size.  However,  when  particle  size  is  smaller  than  the 
wavelength  of  emitted  radiation,  the  particles  radiate  from  their  volume  and  a  size  dependent  emissivity 
could  be  introduced  but  it  is  more  appropriate  to  work  with  a  Rayleigh  low  frequency  expression  for 
particle  absorption/emission  cross  section  at  the  emission  wavelengths  corresponding  to  the  temperature 
elevation.  We  will  consider  only  flare  particle  sizes  greater  than  the  largest  emitted  wavelengths 
~10micrometers  and  so  that  emissivity  does  not  depend  on  particle  size. 

We  turn  our  attention  to  the  ratio  of  conductive  versus  radiative  heat  losses  as  a  function 
of  particle  size,  emissivity  and  temperature.  Conductive  heat  losses  are  4n^K  ATf  and  radiative  heat  losses 
are  ax'*  so  that  the  tatio  of  radiative  over  conductive  heat  loses  at  any  specified  particle 

temperature  is  proportional  to  particle  size.  Setting  this  ratio  equal  to  one  we  find  the  particle  size  feq 
producing  equal  conductive  and  radiative  heat  losses, 

4^  =  l^^T|(s<yr) 

where  T  =  Tambient+AT.  At  temperature  elevations  small  compared  to  ambient  temperature  this  particle 
size  increases  in  proportion  to  AT.  For  example  at  an  emissivity  equal  to  one  and  at  a  3°K  temperature 
elevation  above  an  ambient  temperature  of  300‘TC  the  particle  size  giving  equal  conductive  and  radiative 
heat  losses  is 

^eq(303°K)  =  (6.05X10'^cal/cm/s/^)(4.19j/cal)(3'TC)/((5.67X10-*^j/s/cmW)(303‘lC )'') 

=  0.0159cm. 

Eventually  a  maximum  size  is  reached  at  AT=Tanibient/3  where 

feq  (max)  =  (6.05X10-^cal/cm/s/’K)(4.19J/cal)(100“K)/((5.67X10-'^J/s/cmW)(400°K)'') 

=  0.175cm 

at  an  ambient  temperature  of  300°K.  Finally  for  large  temperature  elevations  this  particle  size  decreases 
in  proportion  to  1/T^,  for  example  at  a  temperature  elevation  of  1500°K, 

^eq(1800"K)  =  (6.05X10-^cal/cm/s/®K)(4.19J/cal)(1500°K)/((  5.67X10-*^J/s/cmW)(1800°K  )^) 

=  0.00639cm 

By  working  with  larger  and  hotter  flare  particles  we  can  increase  the  portion  of  heat 
going  into  cloud  emission  rather  than  heat  conduction  thereby  reducing  contrast  and  signal  to  noise 
transmittance.  By  increasing  the  1800®K  flare  particle  size  in  the  last  expression  two  orders  of  magnitude 
to  0.639cm  diameter  we  can  reduce  conductive  losses  to  1%  of  the  radiative  losses.  But  of  course  such  a 
large  sphere  will  not  stay  airborne  for  long  and  heat  convective  losses  would  also  have  to  be  considered. 
Gravitational  settling  velocities  of  large  particles  can  be  reduced  by  shaping  the  particles  into  tabular 
forms  such  as  disks  or  flakes.  Particle  shapes  having  the  same  surface  area  will  lose  approximately  the 
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same  heat  and  reach  the  same  equilibrium  temperature  as  long  as  particle  heating  and  emission  are  surface 
area  dependent  and  not  volume  dependent.  For  example  thin  disks  roughly  about  0.9cm  in  diameter  can 
be  expected  to  perform  like  the  0.639cm  diameter  sphere  mentioned  but  with  an  adjustable  gravitational 
settling  velocity  and  an  adjustable  surface  chemical  reaction  lifetime  depending  on  the  disk  thickness 
chosen. 


9.  FLARE  PARTICLE  GRAVITATIONAL  SETTLING  VELICITY 

The  settling  velocity  of  disk  shaped  flare  particles  can  be  found  by  looking  at  their  drag 
coefficients  in  the  laminar,  transition  and  turbulent  flow  regimes.  The  gravitational  settling  velocity  is 
then  computed  using  the  equation  defining  the  drag  coefficient 


V.  -25S- 
ic^Ap 

where  m  is  the  mass  of  the  flare  particle  and  A  is  the  cross  sectional  area  of  the  particle  transverse  to  air 
flow.  The  disk  mass  is 


m  =  pr—dh 

for  a  disk  of  density  Pf,  diameter  d  and  thickness  t.  In  the  laminar  flow  regime  the  drag  coefficients  are 

for  motion  parallel  to  the  disk  symmetry  axis  and  — -  for  motion  along  any  equatorial  axis® 
;rRe  3;TtRe 

(perpendicular  to  the  symmetry  axis)  where  the  Reynolds  number  is  defined  in  terms  of  the  major 
dimension 


In  the  turbulent  flow  regime'®’”  the  drag  coefficient  approaches  a  value  of  1  and  the  orientation 
distribution  moves  from  a  glide-tumbling  motion  ultimately  to  complete  random'^  at  very  large  Reynolds 
numbers.  A  simple  approximation  throughout  the  transition  flow  regime  is  a  simple  combination  of  these 
two  solutions 


;rRe 


3;rtRe 


for  V  II  symmetry  axis 


for  V  L  symmetry  axis. 


The  settling  velocities  in  the  laminar  flow  regime  are 


V  =  ^  for  V  II  symmetry 

48  7 


axis 
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V  =  for  V  ±  symmetry  axis 

32tj 


where  tj  =  1  .SSxlO"^  g/cm/s  is  the  viscosity  of  air.  In  the  turbulent  flow  regime  the  gravitational  settling 
velocities  are 


v  = 


V  3p 


for  V  II  symmetry  axis 


V  = 


for  V 1  symmetry  axis. 


Looking  at  the  0.9  cm  diameter  flare  disk  mentioned  earlier  which  puts  most  of  its  heat 
loss  into  radiation  rather  than  conduction,  at  a  thickness  of  20  micrometers  and  a  density  of  2  g/cm^  it 
would  have  a  glide  tumbling  settling  velocity  of  approximately 


Utp^g  ^  1^20x1 0~^cmX2g/cm^)(980cm^ 
3p  ~  ]]  3(0.001 3g/cm^) 

=  63.4  cm/s 


corresponding  to  a  Reynolds  number 


dv  (0.9cm)(63.4cni/s) 

ivC  ■— - — - r - —  3  /  J 

p  0.1 52cm  vs 

which  is  consistent  with  the  chosen  drag  coefficient  and  orientation. 


1 0.  FLARE  PARTICLES  INCORPORATED  INTO  AN  ATTENUATING  SCATTERING 

AEROSOL 

When  considering  flare  particles  in  smoke,  again  the  four  flux  theory  can  be  used  to 
relate  smoke  cloud  contrast  transmittance  and  signal  to  noise  transmittance  to  aerosol  extinction 
coefficient,  single  scatter  albedo,  asymmetry  factor,  concentration  path  length  product  and  cloud 
temperature.  If  we  consider  the  flare  particles  to  be  absorbing  and  therefore  emitting  while  the  smoke 
particles  are  purely  scattering  then  the  temperature  of  the  flare  particles  can  be  adjusted  independent  of 
the  smoke  particle  temperature  and  air  temperature  which  can  remain  at  ambient  air  temperatures  if  the 
flare  particle  size  is  large  enough  to  lose  heat  predominantly  by  radiation  rather  than  conduction  into  the 
ambient  air.  Also  it  is  a  simplification  to  consider  the  smoke  particles  as  purely  scattering  because  they 
will  not  be  heated  by  the  radiation  given  off  by  the  flare  particles.  The  cloud  brightness  will  depend  on 
the  concentration  of  flare  particles. 

Radiation  from  the  cloud  surface  is 


AH 


cloud 
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where  Sc  is  cloud  surface  area,  Tp  is  flare  temperature  and  Scjoud  =  1  -  R  -  Tj  -  Tu,  is  cloud  emissivity.  In 
the  limit  of  an  opaque  cloud  the  emissivity  becomes 


w=i-R  =  i - ^ - 

f(i-g) 


This  is  less  than  the  radiation  that  would  result  if  flare  particles  did  not  radiatively  interact  by  absorbing 
radiation  from  other  flare  particles  and  self  absorption  of  scattered  radiation. 

The  single  scatter  albedo  of  an  aerosol  cloud  containing  absorbing/emitting  flare  particles  embedded 
within  a  purely  scattering  aerosol  which  does  not  heat  up  directly  due  to  radiation  is  found  by  looking  at 
the  fraction  of  radiation  absorbed  (1  -  a)  on  average  by  the  mix  of  flare  and  scattering  particles 


\-a 


CFQ!p(l-Op)  +  C,a,(l-a,) 

CpOp  + 


where  Cp,  ap  and  op  is  the  concentration,  extinction  coefficient,  and  single  scatter  albedo  of  the  flare 
particles  and  Cs,  cts  and  Os  are  the  same  for  the  scattering  particles.  We  note  that  =  1  and  because  the 
flare  particles  are  much  larger  than  the  scattering  smoke  particles  CpiS^p  «:  so  that 


The  absorptivity/emissivity  of  a  black  surface  is  unity  permitting  a  flare  particle  to  emit  the  greatest 
power.  For  such  a  black  particle  there  will  be  no  scatter  cross  section  component  due  to  reflection  or 
refraction  but  there  will  be  a  scatter  component  due  to  diffraction  equal  to  the  absorption  cross  section  so 
we  set  Op  =  1/2  to  maximize  flare  output.  For  flake  or  disk  shaped  flare  and  scattering  particles 


a(mVcm^) 


1 

t  (micrometers) 


where  t  is  the  flare  particle  thickness  tp  or  scattering  smoke  particle  thickness  ts  in  micrometers  and  the 
cloud  average  single  scatter  albedo  becomes 


0  =  1- 


2C,tp 


The  parameter  F  used  in  some  of  the  plots  is  defined  as  the  fraction  of  radiation  absorbed/emitted  on 
average 


F  =  l-o  = 


Cpt, 

2C,tp 


and  therefore  o  =  1  -  F.  Because  each  flare  particle  absorbs  radiation  originating  from  other  flare 
particles  as  well  as  self  emitted  radiation  scattered  back  in  its  direction,  the  power  required  to  maintain 
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flare  temperature  at  any  given  level  is  not  as  great  as  it  would  be  if  the  flare  particles  were  far  apart  and 
acting  independently.  The  surface  area  of  a  spherical  cloud  is  where  D  is  the  cloud  diameter  and  the 
combined  surface  area  S?  of  all  flare  particles  is  the  cloud  volume  ;zDV6  multiplied  by  the  flare  particle 


concentration  Cp  multiplied  by  the  flare  surface  area  over  volume  ratio  2 


r^rdn 

/ 

4  J 

/ 

1  4  J 

tp  -  2/tp 


The  value  of  Cp  can  be  found  from  the  value  of  F  chosen  in  combination  with  the  flare  temperature  and 
cloud  optical  depth  from  the  plots  giving  the  desired  reduction  in  contrast  and  signal  to  noise 
transmittances. 


Cp  - 


p 

ts 


where  the  optical  depth  of  the  cloud  is 


h^iC^a^+C^a^yL 


because  »  CpOi: 


Substituting  this  value  into  the  equation  relating  Cp  to  Cs  and  F 


and  the  total  flare  surface  area  becomes 


Sf  =  |?d^ 


/  o  ^ 


2bt„  F 


Letting  L  =  D 


2;r 


bF 


and  the  radiation  emitted  by  the  total  number  of  flare  particles  acting  independently  is 


18 


The  ratio  f  of  radiation  emitted  by  a  cloud  of  high  emissivity  zjiare «  1  flare  particles  embedded 

within  a  scattering  cloud  divided  by  the  radiation  emitted  by  those  flare  particles  dispersed,  so 

that  they  do  not  interact  or  absorb  their  own  emitted  radiation  due  to  scatter,  is 

j.  (1-R-T.-T.)>rtf 

(l)^D’  b  F 

3(I-R-T.-T.) 

2  bF 

If  a  disk  shaped  cloud  had  been  chosen  y  would  be  less  than  or  equal  to  one  as  it  should  be.  Significant 
reductions  below  this  value  would  occur  only  for  F  >  0.01  and  b  >  3.  The  flare  power  required  per  unit 
volume  of  flare  particles  is  found  by  taking  the  radiation  emitted  by  the  cloud  and  dividing  it  by  the 
volume  of  flare  particles  in  the  cloud 

AH,,^  (l-R-Ts-T,);rDVTp^  BQ-R-Ts -TjtxTp^ 


^-dOcp 

f2btpF^ 

U  ) 

U  J 

D  J 

which  is  on  the  order  of  800  W/cm^  at  Tp  =  1000°K  and  tp  =  0.01cm  over  the  parameter  range  of  interest 
that  significantly  reduces  contrast  and  signal  to  noise  transmittance;  0.001<F<0.1  and  l<b:<3. 


1 1 .  DISCUSSION  OF  FIGURES 

Figures  1-46  address  warm  clouds  where  the  air  temperature  is  close  to  the  aerosol 
particle  temperature.  Contrast  transmittance  is  plotted  in  Figures  1-10.  At  several  values  of  single  scatter 
albedo  and  target  over  background  temperature  ratio  it  appears  as  a  function  of  smoke  optical  depth  and 
cloud  over  background  temperature  ratio  in  Figures  1-4.  At  several  values  of  optical  depth  and  target 
over  background  temperature  ratio  it  appears  as  a  function  of  single  scatter  albedo  and  cloud  over 
background  temperature  ratio  in  Figures  5-8.  All  remaining  plots  assume  that  target  and  background 
temperatures  differ  by  only  a  few  degrees.  Contrast  transmittance  at  two  cloud  over  background 
temperature  ratios  appears  as  a  function  of  smoke  optical  depth  and  single  scatter  albedo  in  Figures  9  and 
10.  Contrast  transmittance  ratio  (contrast  transmittance  at  the  elevated  cloud  temperature  divided  by 
contrast  transmittance  at  ambient  cloud  temperature)  is  computed  in  Figures  11-17.  It  is  computed  in 
Figures  1 1-13  as  a  function  of  cloud  over  background  temperature  ratio  and  optical  depth  at  three  single 
scatter  albedos.  It  is  computed  in  Figures  14  and  15  as  a  Action  of  cloud  over  background  temperature 
ratio  and  single  scatter  albedo  at  two  optical  depths.  It  is  computed  in  Figures  16  and  17  as  a  function  of 
optical  depth  and  single  scatter  albedo  at  two  cloud  over  background  temperature  ratios.  The  ratio  of  the 
derivative  of  contrast  transmittance  with  respect  to  cloud  temperature  over  the  derivative  of  contrast 
transmittance  with  respect  to  optical  depth  is  shown  in  Figures  1 8-24.  This  ratio  is  equal  to  the  ratio  of 
the  optical  depth  increment  that  reduces  contrast  transmittance  by  the  same  amount  that  a  1  degree  Kelvin 
rise  in  cloud  temperature  reduces  contrast  transmittance.  This  optical  depth  increment  is  plotted  as  a 
function  of  optical  depth  and  single  scatter  albedo  in  Figures  18  and  19,  as  a  function  of  cloud 
temperature  and  single  scatter  albedo  in  Figures  20  and  21  and  as  a  function  of  cloud  temperature  and 
optical  depth  in  Figures  22-24.  Signal  to  noise  transmittance  is  shown  in  Figures  25-3 1 .  It  is  computed 
as  a  function  of  optical  depth  and  cloud  over  background  temperature  ratio  in  Figures  25-27,  as  a  function 
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of  single  scatter  albedo  and  cloud  over  background  temperature  ratio  in  Figures  28  and  29  and  as  a 
function  of  optical  depth  and  single  scatter  albedo  in  Figures  30  and  31.  Signal  to  noise  transmittance 
ratio  (signal  to  noise  transmittance  at  the  elevated  cloud  temperature  divided  by  signal  to  noise 
transmittance  at  ambient  cloud  temperature)  is  computed  in  Figures  32-38.  It  is  computed  in  Figures  32- 
34  as  a  function  of  cloud  over  background  temperature  ratio  and  optical  depth  at  three  single  scatter 
albedos.  It  is  computed  in  Figures  35  and  36  as  a  function  of  cloud  over  background  tecnperature  ratio 
and  single  scatter  albedo  at  two  optical  dqjths.  It  is  computed  in  Figures  37  and  38  as  a  ftinction  of 
optical  depth  and  single  scatter  albedo  at  two  cloud  over  backgroxmd  tenperature  ratios.  The  ratio  of  the 
derivative  of  signal  to  noise  transmittance  with  respect  to  cloud  temperature  over  the  derivative  of  signal 
to  noise  transmittance  with  respect  to  optical  depth  is  shown  in  Figures  39-46.  This  ratio  is  equal  to  the 
ratio  of  the  optical  depth  increment  that  reduces  signal  to  noise  transmittance  by  the  same  amoxmt  that  a  1 
degree  Kelvin  rise  in  cloud  temperature  reduces  signal  to  noise  transmittance.  This  optical  depth 
increment  is  plotted  as  a  function  of  optical  depth  and  single  scatter  albedo  in  Figures  39  and  40,  as  a 
function  of  cloud  temperature  and  single  scatter  albedo  in  Figures  41-43  and  as  a  function  of  cloud 
temperature  and  optical  depth  in  Figures  44-46. 

Figures  47-70  address  flare  particles  dispersed  within  a  scattering  smoke  that  is  not 
heated  by  the  flare  particles.  Figures  47-62  show  contrast  transmittance  and  signal  to  noise  transmittance 
as  a  function  of  two  of  the  following  independent  variables  while  holding  the  remaining  variable  fixed ; 
optical  depth,  logarithm  base  10  of  the  ratio  of  flare  over  background  temperature  and  logarithm  base  10 
of  F  or  equivalently  rfi^  which  is  equal  one  minus  the  single  scatter  albedo.  Figures  63-66  show  the 
ratio  of  the  derivative  of  signal  to  noise  transmittance  with  respect  to  the  ratio  of  flare  over  background 
temperature  divided  by  the  derivative  of  signal  to  noise  transmittance  with  respect  to  optical  depth.  This 
is  equal  to  the  optical  depth  increment  required  to  reduce  signal  to  noise  transmittance  by  the  same 
amount  that  a  unit  increment  in  the  ratio  of  flare  over  background  temperature  reduces  signal  to  noise 
transmittance.  Figures  67-70  show  the  ratio  of  the  derivative  of  signal  to  noise  transmittance  with  respect 
to  F  divided  by  the  derivative  of  signal  to  noise  transmittance  with  respect  to  optical  depth.  This  is  equal 
to  the  optical  depth  increment  required  to  reduce  signal  to  noise  transmittance  by  the  same  amoimt  that  a 
unit  increment  in  F  reduces  signal  to  noise  transmittance. 


12.  CONCLUSIONS 

Smoke  optical  depth  can  be  decreased  modestly,  on  the  order  of  0.005  per  degree  Kelvin, 
by  increasing  cloud  temperature  and  emission.  This  number  is  based  on  reductions  in  contrast 
transmittance  and  signal  to  noise  transmittance  over  a  wide  range  of  optical  depths  and  single  scatter 
albedos.  However  the  energy  required  to  raise  the  temperature  of  the  cloud  is  considerable,  about  10,000 
joules  per  gram  of  aerosol  per  degree  Kelvin  at  a  concentration  of  0.13g/m^,  because  of  the  huge  ratio  of 
air  mass  (air  density  is  1300g/m^)  to  aerosol  mass.  In  addition  cloud  buoyancy  quickly  moves  the  cloud 
up  out  of  the  line  of  si^t.  A  ten  meter  diameter  cloud  that  is  3°K  warmer  than  the  surrounding  air  will 
rise  at  a  velocity  of  3.6m/s.  Such  velocities  of  rise  are  in  the  turbulent  flow  regime  and  will  cause  the 
cloud  to  break  up  even  if  it  is  rising  through  relatively  still  air. 

An  alternative  method  to  increase  cloud  emission  has  been  proposed  here  that  is  much 
more  effective.  Cloud  emission  can  be  significantly  increased  by  incorporating  large  tabular  (flake  or 
disk)  shaped  flare  particles  into  a  scattering  smoke  because  a  much  larger  fraction  of  heat  loss  is  put  into 
radiation  rather  than  conduction  into  air  that  does  not  emit  in  atmospheric  windows.  Flare  particle  size, 
shape,  temperature,  concentration  and  the  ratio  of  flare  particles  to  scattering  particles  have  been 
optimized  to  minimize  conductive  heat  loss  and  settling  velocity  while  minimizing  smoke  contrast 
transmittance  and  signal  to  noise  transmittance.  Optical  depth  reductions  on  the  order  of  0.3  for  a  300‘K 
increase  in  flare  particle  temperature  are  possible  over  a  wide  range  of  optical  depths  at  flare  temperatures 
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greater  than  900“K  and  F  values  greater  than  0.00 1 .  We  define  F  to  be  half  the  ratio  of  flare  particle 
concentration  ©ctinction  coefficient  product  divided  by  scattering  particle  concentration  extinction 
coefficient  product.  Optical  depth  reductions  on  the  order  of  0.05  for  a  0.001  increase  in  F  have  been 
computed  over  a  wide  range  of  optical  depths  at  flare  temperatures  above  900°K  and  F  values  greater  than 
0.001 .  Above  flare  particle  temperatures  of 400“K  the  ratio  of  radiation  lose  over  conduction  lose 
increases  with  increasing  particle  size  and  temperature.  For  exanqjle  a  0.9cm  diametCT  flare  disk  at  a 
temperature  of  1800“K  will  lose  99%  of  its  heat  by  radiation  and  have  a  settling  velocity  of  63.4cm/s  at  a 
thickness  of  20  micrometers  and  a  density  of  2g/cc.  To  further  reduce  settling  velocity  flare  thickness  and 
density  might  be  reduced  since  settling  velocity  is  proportional  to  the  square  root  of  the  product  of 
thickness  and  density.  It  is  also  possible  to  create  flare  particle  shapes  that  spin  as  they  fall  to  reduce 
settling  velocity.  Power  production  requirements  for  flare  particles  of  a  given  thickness  are  obtained  by 
looking  at  the  appropriate  figure  and  then  choosing  a  combmation  of  optical  depth,  F  and  flare 
temperature  that  reduces  contrast  transmittance  or  signal  to  noise  transmittance  to  the  desired  level.  For 
example  flare  particles  produce  800w/cc  at  a  ten5)erature  of  1000"K  and  thickness  of  100  micrometers 
over  the  range  of  optical  depths  and  F  values  of  interest. 

In  closing,  certain  limitations  should  be  mentioned.  This  study  has  represented  emitted 
power  as  proportional  to  ten:5)erature  raised  to  the  fourth  power  without  addressing  wavelength 
dependence  which  manifests  itself  as  peak  power  output  moving  progressively  to  shorter  wavelengths  as 
flare  temperature  is  increased.  Also  all  flare  particle  calculations  of  contrast  transmittance  and  signal  to 
noise  transmittance  have  assumed  a  scattering  nonabsorbing  smoke.  Real  smoke  particles  also  absorb  and 
reemit  radiation  in  at  least  portions  of  the  infrared  but  this  effect  could  not  be  studied  using  the  simplified 
radiative  transfer  equations  presented  here.  Finally,  clouds  were  assumed  to  be  homogeneous  walls  or 
blankets  while  in  face  the  clouds  are  inhomogeneoris  and  finite  in  all  dimensions.  This  means  that  effects 
due  to  flux  transported  through  the  sides  of  clouds  were  not  investigated  nor  were  the  effects  of  clutter 
due  to  aerosol  concentration  inhomogeneities  studied. 
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SOURCE  CODE 


Cloud  Infrared  Emission  versus  Attenuation  Codes 

Program  Darpacrl 
g  =  0.75; 

temptargovrbkg  =1.01; 
p  =  a(l  +  g)/2; 
q  =  a(l  -  g)/2; 
z  =  Sqrt[(l  -a)(l-ag)]; 

Tu  =  Exp[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Exp[z  b]  -  Exp[-z  b])/((l  -  p  +  z)E)q)[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

(♦contrast  ratio  for  cloud  at  ten:5)cloudovrbkg  divided  by  contrast  for  cloud  \ 
of  same  properties  at  backgroimd  teii5)erature*) 

CR  =  (Tu(l  +  ten^targovrbkgM)  +  2(R  +  T  -  Tu)  + 

2(1  -  R  -  T))/(Tu(l  +  ten:5)targovrbkg^4)  +  2(R  +  T  -  Tu)  + 

2(1  -  R  -  T)ten:q)cloudovrbkg'^4  ); 
a  =  0.5; 

Plot3D[CR,  {tempcloudovrbkg,  1,  1.5},  {b,  1,  3}, 

PlotLabel  ->  ”CR=CT(Tc)/CT(Tc=Tb)  Tt~Tb,a=0.5", 

AxesLabel  ->  {"Tc/Tb",  "OD",  "CR  "},  PlotPoints  ->  25,  PlotRange  ->  (0, 1}]; 
a  =  0.1; 

Plot3D[CR,  {tenqjcloudovrbkg,  1,  1.5},  {b,  1,  3}, 

PlotLabel  ->  "CR=CT(Tc)/CT(Tc=Tb)  Tt~Tb,a=0.1", 

AxesLabel  ->  {"Tc/Tb",  "OD",  "CR  "},  PlotPoints  ->  25,  PlotRange  ->  {0,  1}]; 
a  =  0.9; 

Plot3D[CR,  {tenqKloudovrbkg,  1,  1.5},  {b,  1,  3}, 

PlotLabel  ->  "CR=CT(Tc)/CT(Tc=Tb)  Tt~Tb,a=0.9", 

AxesLabel  ->  {"Tc/Tb",  "OD",  "CR  "},  PlotPoints  ->  25,  PlotRange  ->  {0,  1}] 
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Program  Darpacrl  1 
g  =  0.75; 

ten5)targovrbkg  =1.01; 
p  =  a(l  +  g)/2; 
q  =  a(l  -  g)/2; 
z  =  Sqrt[(l  -a)(l  -ag)]; 

Tu  =  Exp[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Ejq)[z  b]  -  E3q)[-z  b])/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 
(♦contrast  ratio  for  cloud  at  tenqjcloudovrbkg  divided  by  contrast  for  cloud  \ 
of  same  properties  at  background  ten:5)erature*) 

CR  =  (Tu(l  +  tenptargovrbkgM)  +  2(R  +  T  -  Tu)  + 

2(1  -  R  -  T))/(Tu(l  +  ten:5)targovrbkg^4)  +  2(R  +  T  -  Tu)  + 

2(1  -  R  -  T)tenq)cloudovrbkgM ); 

SNTR  =  CR^5; 
a=  0.5; 

Plot3D[SNTR,  {tempcloudovrbkg,  1, 1.5},  {b,  1,  3}, 

PlotLabel  ->  ”SNTR=SNT(Tc)/SNT(Tc=Tb)  Tt~Tb,a=0.5", 

AxesLabel  ->  {"Tc/Tb",  "OD",  "SNTR  PlotPoints  ->  25, 

PlotRange  ->  {0, 1}]; 
a  =  0.1; 

Plot3D[SNTR,  {ten^loudovrbkg,  1,  1.5},  {b,  1,  3}, 

PlotLabel  ->  "SNTR=SNT(Tc)/SNT(Tc=Tb)  Tt~Tb,a=0.1", 

AxesLabel  ->  {"Tc/Tb",  "OD",  "SNTR  "},  PlotPoints  ->  25, 

PlotRange  ->  {0, 1}]; 
a  =0.9; 

Plot3D[SNTR,  {tempcloudovrbkg,  1, 1.5},  {b,  1,  3}, 

PlotLabel  ->  "SNTR=SNT(Tc)/SNT(Tc=Tb)  Tt~Tb,a=0.9", 

AxesLabel  ->  {"Tc/Tb",  "OD",  "SNTR  "},  PlotPoints  ->  25, 

PlotRange  ->  {0, 1}] 
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Program  Darpacr2 
g  =  0.75; 

tenqjtargovrbkg  =1.01; 
p  =  a(l  +  g)/2; 
q  =  a(l-g)/2; 
z=SqTt[(l  -a)(l  -ag)]; 

Tu  =  Exp[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Exp[z  b]  -  Exp[-z  b])/((l  -  p  +  z)E3q)[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 
(♦contrast  ratio  for  cloud  at  ten5)cloudovrbkg  divided  by  contrast  for  cloud  \ 
of  same  properties  at  background  ten:5)erature*) 

CR  =  (Tu(l  +  ten5)targovrbkg^4)  +  2(R  +  T  -  Tu)  + 

2(1  -  R  -  T))/(Tu(l  +  ten5)targovrbkgM)  +  2(R  +  T  -  Tu)  + 

2(1  -  R  -  T)ten:5)cloudovrbkg'^4 ); 

b=l; 

Plot3D[CR,  {ten:5)cloudovrbkg,  1, 1.5},  {a,  0.0001,  .9999}, 

PlotLabel  ->  "CR=CT(Tc)/CT(Tc=Tb)  Tt~Tb,ly=l", 

AxesLabeL->  {"Tc/Tb",  "a",  "CR "},  PlotPoints  ->  25,  PlotRange  ->  {0, 1}]; 
b=1.5; 

Plot3D[CR,  {teirqxiloudovrbkg,  1,  1.5},  {a,  0.0001,  .9999}, 

PlotLabel  ->  "CR=CT(Tc)/CT(Tc=Tb)  Tt~Tb,b=1.5", 

AxesLabel  ->  {"Tc/Tb",  "a",  "CR "},  PlotPoints  ->  25,  PlotRange  ->  {0,  1}] 
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Program  Darpacr22 
g  =  0.75; 

tenqjtargovrbkg  =  1.01; 
p  =  a(l  +g)/2; 
q  =  a(l  -g)/2; 
z  =  Sqrt[(l  -  a)(l  -  a  g)]; 

Tu  =  Ejq)[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Exp[z  b]  -  Exp[-z  b])/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Ejq3[-z  b]); 

(♦contrast  ratio  for  cloud  at  tenqjcloudovrbkg  divided  by  contrast  for  cloud  \ 
of  same  properties  at  backgroimd  tenq)erature*) 

CR  =  (Tu(l  +  ten:q)targovrbkgM)  +  2(R  +  T  -  Tu)  + 

2(1  -  R  -  T))/(Tu(l  +  tenqjtargovrbkgM)  +  2(R  +  T  -  Tu)  + 

2(1  -  R  -  T)ten5)cloudovrbkgM ); 

SNTR  =  CR^5; 
b  =  l; 

Plot3D[SNTR,  {tenpcloudovrbkg,  1, 1.5},  {a,  .0001,  .9999}, 

PlotLabel  ->  "SNTR=SNT(Tc)/SNT(Tc=Tb)  Tt~Tb,b=l", 

AxesLabel  ->  {"Tc/Tb",  "a",  "SNTR  "},  PlotPoints  ->  25,  PlotRange  ->  {0, 1}]; 
b=1.5; 

Plot3D[SNTR,  {ten^jcloudovrbkg,  1, 1.5},  {a,  .0001,  .9999}, 

PlotLabel  ->  "SNTR=SNT(Tc)/SNT(Tc=Tb)  Tt~Tb,b=1.5", 

AxesLabel  ->  {"Tc/Tb",  "a",  "SNTR  "},  PlotPoints  ->  25,  PlotRange  ->  {0, 1}] 
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Program  DarpacrS 
g  =  0.75; 

ten:5)targovrbkg  =  1.01; 
p  =  a(l  +g)/2; 
q  =  a(l  -g)/2; 
z  =  Sqrt[(l  -a)(l  -ag)]; 

Tu  =  Exp[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Ejq)[z  b]  -  Exp[-z  b])/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 
(♦contrast  ratio  for  cloud  at  ten^cloudovrbkg  divided  by  contrast  for  cloud  \ 
of  same  properties  at  backgroimd  ten^erature*) 

CR  =  (Tu(l  +  ten:5)targovrbkg'^4)  +  2(R  +  T  -  Tu)  + 

2(1  -  R  -  T))/(Tu(l  +  temptargovrbkg^4)  +  2(R  +  T  -  Tu)  + 

2(1  -  R  -  T)ten:5)cloudovrbkgM  ); 
ten^cloudovrbkg  =1.1; 

Plot3D[CR,  {b,  1,  3},  {a,  0.0001,  .9999}, 

PlotLabel  ->  "CR=CT(Tc)/CT(Tc=Tb)  Tt~Tb,Tc/Tb=l.l”, 

AxesLabel  ->  {"b",  "a",  "CR  "},  PlotPoints  ->  25,  PlotRange  ->  {0, 1}]; 
ten5)cloudovrbkg  =  1.5; 

Plot3D[CR,  {b,  1,  3},  {a,  0.0001,  .9999}, 

PlotLabel  ->  "CR=CT(Tc)/CT(Tc=Tb)  Tt~Tb,Tc/Tb=1.5", 

AxesLabel  ->  {"b",  "a",  "CR  "},  PlotPoints  ->  25,  PlotRange  ->  {0,  1}] 
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Program  Darpacr33 
g  =  0.75; 

tenqjtargovrbkg  =1.01; 
p  =  a(l  +  g)/2; 
q  =  a(l  -  g)/2; 
z  =  Sqrt[(l  -a)(l  -ag)]; 

Tu  =  Exp[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Exp[z  b]  -  Exp[-z  b])/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 
(♦contrast  ratio  for  cloud  at  tempcloudovrbkg  divided  by  contrast  for  cloud  \ 
of  same  properties  at  background  tenqjerature*) 

CR  =  (Tu(l  +  tenptargovrbkgM)  +  2(R  +  T  -  Tu)  + 

2(1  -  R  -  T))/(Tu(l  +  temptargovrbkg^4)  +  2(R  +  T  -  Tu)  + 

2(1  -  R  -  T)terapcloudovrbkgM ); 

SNTR  =  CR^5; 
tempcloudovrbkg  =  1.1; 

Plot3D[SNTR,  {b,  1,  3},  {a,  .0001,  .9999}, 

PlotLabel  ->  "SNTR=SNT(Tc)/SNT(Tc=Tb)  Tt~Tb,Tc/Tb=l.l", 
AxesLabel  ->  {"b",  "a",  "SNTR  PlotPoints  ->25,  PlotRange  ->  {0, 1}]; 
tenq)cloudovrbkg  =  1.5; 

Plot3D[SNTR,  {b,  1,  3},  {a,  .0001,  .9999}, 

PlotLabel  ->  "SNTR=SNT(Tc)/SNT(Tc=Tb)  Tt~Tb,Tc/Tb=1.5", 
AxesLabel  ->  {"b",  "a",  "SNTR  "},  PlotPoints  ->  25,  PlotRange  ->  {0, 1}] 
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Program  Darpactl 
g  =  0.75; 
p  =  a(l  +  g)/2; 
q  =  a(l-g)/2; 
z  =  Sqrt[(l  -  a)(l  -  a  g)]; 

Tu  =  Ejq)[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(E3q)[z  b]  -  Exp[-z  b])/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

(♦Note  that  contrast  transmittance  is  just  contrast  ratio  - 
contrast  with  cloud/without  cloud*) 

CT  =  Tu(ten5)targovrbkgM  +  l)/(Tu(ten:5)targovrbkg^4  +  1)  +  2(R  +  T  -  Tu)  + 
2(1  -  R  -  T)ten:5)cloudovrbkgM ); 

a  =  0.5; 

ten:5)targovrbkg  =  1.0001; 

gl  =Plot3D[CT,  {b,  1, 3},  {terapcloudovrbkg,  1, 1.5}, 

PlotLabel  ->  "Contrast  Transmittance  Tt~Tb,a=0.5", 

AxesLabel  ->  {"OD",  "Tc/Tb",  "CT  "},  PlotPoints  ->  25]; 
a  =  0.1; 

ten5)targovrbkg  =  1.0001; 

gl  =Plot3D[CT,  {b,  1, 3},  {tempcloudovrbkg,  1, 1.5}, 

PlotLabel  ->  "Contrast  Transmittance  Tt~Tb,a=0.1", 

AxesLabel  ->  {"OD",  "Tc/Tb",  "CT  "},  PlotPoints  ->  25]; 
a  =  0.9; 

ten:5)targovrbkg  =  1.0001; 

gl  =Plot3D[CT,  {b,  1,  3},  {tenqjcloudovrbkg,  1, 1.5}, 

PlotLabel  ->  "Contrast  Transmittance  Tt~Tb,a=0.9", 

AxesLabel  ->  ("OD",  "Tc/Tb",  "CT  "},  PlotPoints  ->  25] 
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Program  Darpactl  1 
g  =  0.75; 
p  =  a(l  +  g)/2; 
q  =  a(l  -  g)/2; 
z  =  Sqrt[(l-a)(l-ag)]; 

Tu  =  E3q)[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Exp[z  b]  -  Exp[-z  b])/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

(*Note  that  contrast  transmittance  is  just  contrast  ratio  - 
contrast  with  cloud/without  cloud*) 

CT  =  Tu(ten]ptargovrbkg'^4  +  l)/(Tu(ten5)targovrbkgM  +  1)  +  2(R  +  T  -  Tu)  + 
2(1  -  R  -  T)tenpcloudovrbkgM  ); 

SNT  =  (TuCT)^5; 
a  =  0.5; 

ten:5)targovrbkg  =  1.0001; 

gl  =Plot3D[SNT,  {b,  1,  3},  {teii^cloudovrbkg,  1, 1.5}, 

PlotLabel  ->  "S/N  Transmittance  Tt~Tb,a=0.5", 

AxesLabel  ->  {"OD",  "Tc/Tb",  "SNT  "},  PlotPoints  ->  25]; 
a  =  0.1; 

tenqjtargovrbkg  =  1.0001; 

gl  =Plot3D[SNT,  {b,  1,  3},  {teirqjcloudovrbkg,  1, 1.5}, 

PlotLabel  ->  "S/N  Transmittance  Tt~Tb,a=0.r', 

AxesLabel  ->  {"OD",  "Tc/Tb",  "SNT  "},  PlotPoints  ->  25]; 
a  =  0.9; 

teiEptargovrbkg  =  1 .000 1 ; 

gl  =Plot3D[SNT,  {b,  1,  3},  {teirpcloudovrbkg,  1, 1.5}, 

PlotLabel  ->  "S/N  Transmittance  Tt~Tb,a=0.9", 

AxesLabel  ->  {"OD",  "Tc/Tb",  "SNT  "},  PlotPoints  ->  25] 
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Program  Darpact2 
g  =  0.75; 
p  =  3(1+  g)/2; 
q  =  a(l-g)/2; 
z  =  Sqrt[(l  -  a)(l  -  a  g)]; 

Tu  =  Exp[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Ejq)[z  b]  -  Exp[-z  b])/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

(*Note  that  contrast  transmittance  is  just  contrast  ratio  - 
contrast  with  cloud/without  cloud*) 

CT  =  Tu(ten:5)targovrbkgM  +  l)/(Tu(ten:^targovrbkgM  +  1)  +  2(R  +  T  -  Tu)  + 
2(1  -  R  -  T)ten:5)cloudovrbkgM  ); 

b=l; 

ten:5)targovrbkg  =  1.0001; 

Plot3D[CT,  {tenqjcloudovrbkg,  1, 1.5},  {a,  .999,  0.0001}, 

PlotL^l  ->  "Contrast  Transmittance  Tt~Tb,OD=l", 

AxesLabel  ->  {"Tc/Tb",  "a",  "CT "},  PlotPoints  ->  25]; 

Show[%,  PlotRmge  ->  {0,  .4}]; 
b=1.5; 

ten5)targovrbkg  =  1.0001; 

Plot3D[CT,  {tempcloudovrbkg,  1, 1.5},  {a,  .999,  0.0001}, 

PlotLabel  ->  "Contrast  Transmittance  Tt~Tb,OD=1.5", 

AxesLabel  ->  {"Tcm)",  "a",  "CT "},  PlotPoints  ->  25]; 

Show[%,  PlotRmge  ->  {0,  .25}] 
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Program  Darpact22 
g  =  0.75; 
p  =  a(l  +  g)/2; 
q  =  a(l  -  g)/2; 
z  =  Sqrt[(l  -a)(l  -ag)]; 

Tu  =  Exp[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Exp[z  b]  -  Ejq)[-z  b])/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

(♦Note  that  contrast  transmittance  is  just  contrast  ratio  - 
contrast  with  cloud/without  cloud*) 

CT  =  Tu(ten:q)targovrbkg^4  +  l)/(Tu(ten5)targovrbkgM  +  1)  +  2(R  +  T  -  Tu)  + 
2(1  -  R  -  T)tenq)cloudovrbkgM ); 

SNT  =  (TuCT)^5; 
b=l; 

temptargovrbkg  =  1.0001; 

gl  =Plot3D[SNT,  {tenq)cloudovrbkg,  1, 1.5},  {a,  .9999,  0.00001}, 

PlotLabel  ->  "S/N  Transmittance  Tt~Tb,OD=l", 

AxesLabel  ->  {"Tc/Tb",  "a",  "SNT  "},  PlotPoints  ->  25]; 

Show[%,  PlotRange  ->  {0,  .4}]; 
b=1.5; 

gl  =Plot3D[SNT,  {ten:5)cloudovrbkg,  1,  1.5},  {a,  .9999,  0.00001}, 

PlotLabel  ->  "S/N  Transmittance  Tt~Tb,OD=1.5", 

AxesLabel  ->  {"Tcm)",  "a",  "SNT  "},  PlotRange  ->  {0,  .25}, 

PlotPoints  ->  25] 
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Program  Darpact3 
g  =  0.75; 
p  =  a(l  +  g)/2; 
q  =  a(l  -  g)/2; 
z  =  Sqrt[(l  -  a)(l  -  a  g)]; 

Tu  =  Exp[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Exp[z  b]  -  Exp[-z  b])/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

(♦Note  that  contrast  transmittance  is  just  contrast  ratio  - 
contrast  with  cloud/without  cloud*) 

CT  =  Tu(ten:5)targovrbkgM  +  l)/(Tu(teinptargovrbkg^4  +  1)  +  2(R  +  T  -  Tu)  + 
2(1  -  R  -  T)ten5)cloudovrbkg'^4  ); 
tenqjcloudovrbkg  =  1.1; 
ten^targovrbkg  =  1.0001; 
gl  =Plot3D[CT,  {b,  1, 3},  {a,  0.00001,  .99999}, 

PlotLabel  ->  "Contrast  Transmittance  Tt~Tb,Tc/Tb=l.r', 

AxesLabel  ->  {"OD",  "a",  "CT  "},  PlotPoints  ->  25]; 
ten:5)cloudovrbkg  =  1.5; 

gl  =Plot3D[CT,  {b,  1, 3},  {a,  0.00001,  .99999}, 

PlotLabel  ->  "Contrast  Transmittance  Tt~Tb,Tc/Tb=1.5", 

AxesLabel  ->  {"OD",  "a",  "CT  "},  PlotPoints  ->  25,  PlotRange  ->  {0,  .35}] 
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Program  Darpact33 
g  =  0.75; 
p  =  a(l  +  g)/2; 
q  =  a(l-g)/2; 
z=Sqrt[(l-a)(l-ag)]; 

Tu  =  Exp[-b]; 

T  =  2z/((l  -  p  +  z)E)q)[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Exp[z  b]  -  Exp[-z  b])/((l  -  p  +  z)Ejq)[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

(*Note  that  contrast  transmittance  is  just  contrast  ratio  - 
contrast  with  cloud/vsathout  cloud*) 

CT  =  Tu(ternptargovrbkg^4  +  l)/(Tu(ten:5)targovrbkgM  +  1)  +  2(R  +  T  -  Tu)  + 
2(1  -  R  -  T)tenq)cloudovrbkgM ); 

SNT  =  (TuCT)^5; 

ten^cloudovrbkg  =1.1; 

temptargovrbkg  =  1.0001; 

gl  =Plot3D[SNT,  {b,  1,  3},  {a,  .9999,  0.00001}, 

PlotLabel  ->  "S/N  Transmittance  Tt~Tb,Tc/TT3=l.l", 

AxesLabel  ->  {"OD",  "a",  "SNT  "},  PlotPoints  ->  25]; 
ten^cloudovrbkg  =1.5; 

gl  =Plot3D[SNT,  {b,  1,  3},  {a,  .9999,  0.00001}, 

PlotLabel  ->  "S/N  Transmittance  Tt~Tb,Tc/TT3=1.5", 

AxesLabel  ->  {"OD",  "a",  "SNT  "},  PlotPoints  ->  25] 
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Program  Darpaemitl 

(♦found  plots  are  independent  of  g  as  expected  for  homogeneous  illumination*) 

g  =  0.75; 

tback  =  300; 

ttarg  =  301; 

p  =  a(l  +  g)/2; 

q  =  a(l  -  g)/2; 

z  =  Sqrt[(l  -a)(l  -ag)]; 

Tu  =  Exp[-b]; 

T  =  2z/((l  -  p  +  z)E3q)[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(E5q)[z  b]  -  Exp[-z  b])/((l  -  p  +  z)E3q)[z  b]  -  (1  -  p  -  z)Ejq)[-z  b]); 

CC  =  Tu(ttarg^4  -  tbackM)/(Tu(ttarg^4  +  tbackM)  +  2  tback^4(R  +  T  -  Tu)  + 
2(l-R-T)tcloudM); 
dd  =  D[CC,  tcloud]/D[CC,  b]; 
tcloud  =  300; 

(*FindMinimum[(CC  -  .02)^2,  {b,  1}]*) 

Plot3D[dd,  {b,  0.0001,  3},  {a,  0.0001,  0.9999}, 

PlotLabel  ->  "(dC/dTc)/(dC/dOD)  Tt~Tb~Tc~300K", 

AxesLabel  ->  {"OD", "  a",  "dOD/dTc  "},  PlotPoints  ->  25, 

PlotRange  ->  {0,  .015}]; 
tcloud  =  400; 

Plot3D[dd,  {b,  0.0001,  3},  {a,  0.0001,  0.9999}, 

PlotLabel  ->  "(dC/dTc)/(dC/dOD)  Tt~Tb~300,To^00", 

AxesLabel  ->  {"OD", "  a",  "dOD/dTc  "},  PlotPoints  ->  25, 

PlotRange  ->  {0,  .015}] 
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Program  Darpaemitl  1 

(♦foxmd  plots  are  independent  of  g  as  expected  for  homogeneous  illumination*) 

g  =  0.75; 

tback  =  300; 

ttarg  =  301; 

p  =  a(l  +g)/2; 

q  =  a(l  -g)/2; 

z  =  Sqrt[(l  -  a)(l  -  a  g)]; 

Tu  =  Ejq)[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Exp[z  b]  -  Exp[-z  b])/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

CC  =  Tu(ttargM  - 

tbackM)/(Tu(ttargM  +  tbackM)  +  2  tbackM(R  +  T  -  Tu)  + 

2(1  -  R  -  T)tcloudM  )^5; 
dd  =  D[CC,  tcloud]/D[CC,  b]; 
tcloud  =  300; 

(*FindMinimum[(CC  -  .02)^2,  {b,  1}]*) 

Plot3D[dd,  (b,  0.0001,  3},  {a,  0.0001,  0.9999}, 

PlotLabel  ->  "(dSN/dTc)/(dSN/dOD)  Tt~Tb~Tc~300K", 

AxesLabel  ->  {"OD",  "  a",  "dOD/dTc  "},  PlotPoints  ->  25, 

PlotRange  ->  {0,  .007}]; 
tcloud  =  600; 

Plot3D[dd,  {b,  0.0001,  3},  {a,  0.0001,  0.9999}, 

PlotLabel  ->  "(dSN/dTc)/(dSN/dOD)  Tt~Tb~300,Tc~600", 

AxesLabel  ->  {"OD",  "  a",  "dOD/dTc  "},  PlotPoints  ->  25, 

PlotRange  ->  {0,  .007}] 
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Program  Darpaemit2 

(*found  plots  are  independent  of  g  as  expected  for  homogeneous  illumination*) 

g  =  0.75; 

tback  =  300; 

ttarg  =  301; 

p  =  a(l  +  g)/2; 

q  =  a(l-g)/2; 

z  =  Sqrt[(l  -  a)(l  -  a  g)]; 

Tu  =  E}q)[-b]; 

T  =  2z/((l  -  p  +  z)Ejq)[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Exp[z  b]  -  Exp[-z  b])/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

CC  =  Tu(ttarg'^4  -  tback'^4)/(Tu(ttargM  +  tbackM)  +  2  tback^4(R  +  T  -  Tu)  + 
2(l-R-T)tcloudM); 
dd  =  D[CC,  tcloud]/D[CC,  b]; 
b=l; 

(*FindMinimum[(CC  -  .02)^2,  {b,  1}]*) 

Plot3D[dd,  {tcloud,  300,  600},  {a,  0.0001,  0.9999}, 

PlotLabel  ->  "(dC/dTc)/(dC/dOD)  OD=l,Tt~Tb~300K", 

AxesLabel  ->  {"Tc",  "  a",  "dOD/dTc  "},  PlotPoints  ->  25, 

PlotRange  ->  {0,  .015}]; 
b  =  2; 

Plot3D[dd,  {tcloud,  300,  600},  {a,  0.0001,  0.9999}, 

PlotLabel  ->  "(dC/dTc)/(dC/dOD)  OD=2,Tt"'Tb~300K", 

AxesLabel  ->  {"Tc",  "  a",  "dOD/dTc  "},  PlotPoints  ->  25, 

PlotRange  ->  {0,  .015}] 


51 


Program  Darpaemit22 

(*foxmd  plots  are  independent  of  g  as  expected  for  homogeneous  dlumination*) 

g  =  0.75; 

tback  =  300; 

ttarg  =  301; 

p  =  a(l  +  g)/2; 

q  =  a(l  -  g)/2; 

z=Sqrt[(l-a)(l-ag)]; 

Tu  =  E3qp[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Exp[z  b]  -  Exp[-z  b])/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

CC  =  Tu(ttargM- 

tback''4)/(Tu(ttargM  +  tback^4)  +  2  tback'^4(R  +  T  -  Tu)  + 
2(l-R-T)tcloudM)^5; 
dd  =  D[CC,  tcloud]/D[CC,  b]; 
b=  1; 

(*FindMinimum[(CC  -  .02)^2,  {b,  1}]*) 

Plot3D[dd,  {tcloud,  300,  600},  {a,  0.0001,  0.9999}, 

PlotLabel  ->  "(dSN/dTc)/(dSN/dOD)  OD=l,Tt~Tb~300K", 

AxesLabel  ->  {"Tc", "  a",  "dOD/dTc  "},  PlotPoints  ->  25, 

PlotRange  ->  {0,  .007}]; 
b  =  2; 

(*FindMmimum[(CC  -  .02)^2,  {b,  1}]*) 

Plot3D[dd,  {tcloud,  300,  600},  {a,  0.0001,  0.9999}, 

PlotLabel  ->  "(dSN/dTc)/(dSN/dOD)  OD=2,Tt~Tb~300K", 

AxesLabel  ->  {"Tc", "  a",  "dOD/dTc  "},  PlotPoints  ->  25, 

PlotRange  ->  {0,  .007}]; 
b  =  3; 

(*FindMinimum[(CC  -  .02)^2,  {b,  1}]*) 

Plot3D[dd,  {tcloud,  300,  600},  {a,  0.0001,  0.9999}, 

PlotLabel  ->  "(dSN/dTc)/(dSN/dOD)  OD=3,Tt~Tb~300K", 

AxesLabel  ->  {"Tc",  "  a",  "dOD/dTc  "},  PlotPoints  ->  25, 

PlotRange  ->  {0,  .007}] 
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Program  DaipaemitS 

(*found  plots  are  independent  of  g  as  expected  for  homogeneous  illumination*) 

g  =  0.75; 

tback  =  300; 

ttarg  =  301; 

p  =  a(l  +g)/2; 

q  =  a(l-g)/2; 

z  =  Sqrt[(l  -a)(l-ag)]; 

Tu  =  E)^[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Exp[z  b]  -  Ejq)[-z  b])/((l  -  p  +  z)E;qp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

CC  =  Tu(ttarg^4  -  tbackM)/(Tu(ttargM  +  tteckM)  +  2  tbackM(R  +  T  -  Tu)  + 
2(l-R-T)tcloudM); 
dd  =  D[CC,  tcloud]/D[CC,  b]; 
a  =  .1; 

(*FindMinimum[(CC  -  .02)'^2,  {b,  1}]*) 

Plot3D[dd,  {tcloud,  300, 600},  {b,  0, 3}, 

PlotLabel  ->  "(dC/dTc)/(dC/dOD)  a=0.1,Tt~lV300K", 

AxesLabel  ->  {"Tc",  "  OD",  "dOD/dTc  ”},  PlotPoints  ->  25, 

PlotRange  ->  {0,  .015}]; 
a  =  .5; 

Plot3D[dd,  {tcloud,  300,  600},  {b,  0, 3}, 

PlotLabel  ->  "(dC/dTc)/(dC/dOD)  a=0.5,Tt~TV'300K", 

AxesLabel  ->  {"Tc",  "  OD",  "dOD/dTc  "},  PlotPoints  ->  25, 

PlotRange  ->  {0,  .015}]; 
a  =  .9; 

Plot3D[dd,  {tcloud,  300,  600},  {b,  0, 3}, 

PlotLabel  ->  "(dC/dTc)/(dC/dOD)  a=0.9,Tt~Tb~300K", 

AxesLabel  ->  {"Tc",  "  OD",  "dOD/dTc  "},  PlotPoints  ->  25, 

PlotRange  ->  {0,  .015}] 
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Program  Darpaemit33 

(♦found  plots  are  independent  of  g  as  expected  for  homogeneous  illumination*) 

g  =  0.75; 

tback  =  300; 

ttarg  =  301; 

p  =  a(l  +  gy2; 

q  =  a(l  -  g)/2; 

z  =  Sqrt[(l  -  a)(l  -  a  g)]; 

Tu  =  Exp[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(E3q)[z  b]  -  Exp[-z  b])/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

CC  =  Tu(ttargM  - 

tbackM)/(Tu(ttarg^4  +  tback^4)  +  2  tback^4(R  +  T  -  Tu)  + 
2(l-R-T)tcloudM)^5; 
dd  =  D[CC,  tcloud]/D[CC,  b]; 
a  =  0.1; 

(*FmdMinimum[(CC  -  .02)'^2,  {b,  1}]*) 

Plot3D[dd,  {tcloud,  300, 600},  {b,  0, 3}, 

PlotLabel  ->  "(dSN/dTc)/(dSN/dOD)  a=0.1,Tt~Tb~300K", 

AxesLabel  ->  {"Tc", "  OD",  "dOD/dTc  "},  PlotPoints  ->  25, 

PlotRange  ->  {0,  .007}]; 
a  =  0.5; 

(*FindMinimum[(CC  -  .02)^2,  {b,  1}]*) 

Plot3D[dd,  {tcloud,  300,  600},  {b,  0,  3}, 

PlotLabel  ->  "(dSN/dTc)/(dSN/dOD)  a=0.5,Tt~Tb~300K", 

AxesLabel  ->  {"Tc",  "  OD",  "dOD/dTc  "},  PlotPoints  ->  25, 

PlotRange  ->  {0,  .007}]; 
a  =  0.9; 

(♦FnidMinimum[(CC  -  .02)^2,  {b,  1}]*) 

Plot3D[dd,  {tcloud,  300, 600},  {b,  0,  3}, 

PlotLabel  ->  "(dSN/dTc)/(dSN/dOD)  a=0.9,Tt~Tb~300K", 

AxesLabel  ->  {"Tc", "  OD",  "dOD/dTc  "},  PlotPoints  ->  25, 

PlotRange  ->  {0,  .007}] 
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Program  Darpaflare  1 

(*ir  flare  emission  within  a  cloud  is  represented  with  a  high  tenqjerature  \ 
high  albedo  cloud  where  1  - 
albedo  is  the  absorption/ 

emission  which  is  just  the  volume  fraction  of  the  flare  /  (flare  + 
white  smoke)  material  multiplied  by  the  ratio  of  flare  absorption  \ 
coefficient  divided  by  smoke  extinction  coefficient  ~ 
radius  smoke  particles  /  Radius  flare  particles*) 
g  =  0.75; 

p  =  a(l  +  g)/2; 
q  =  a(l-g)/2; 
z  =  Sqrt[(l  -a)(l  -ag)]; 

Tu  =  Exp[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Exp[z  b]  -  Exp[-z  b])/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

(*Note  that  contrast  transmittance  is  just  contrast  ratio  - 
contrast  with  cloud/without  cloud*) 

CT  =  Tu(ten:5)targovrbkgM  +  l)/(Tu(ten5)targovrbkg^4  +  1)  +  2(R  +  T  -  Tu)  + 
2(1  -  R  -  TjtenqjcloudovrbkgM ); 

b=  1; 

a=l  -lO^ffi 

tenqjtargovrbkg  =  1,0001; 

Plot3D[CT,  {ten:5)cloudovrbkg,  2,  6},  {ffi  -2,  -5}, 

PlotLabel  ->  "Flare  Contrast  Transmittance  Tt~Tb,OI>=l", 

AxesLabel  ->  {"TfTb",  "log  rfrR",  "CT  PlotPoints  ->  25]; 
b=1.5; 

Plot3D[CT,  {terr5)cloudovrbkg,  2,  6},  {ffi  -2,  -5}, 

PlotLabel  ->  "Flare  Contrast  Transmittance  Tt~Tb,OI>=1.5", 

AxesLabel  ->  {"TfTb",  "log  rfR",  "CT  PlotPoints  ->  25]; 
b  =  2; 

Plot3D[CT,  {ten5)cloudovrbkg,  2, 6},  {ffi  -2,  -5}, 

PlotLabel  ->  "Flare  Contrast  Transmittance  Tt~Tb,OD=2", 

AxesLabel  ->  {"TfrTb",  "log  rfrR",  "CT "},  PlotPoints  ->  25]; 
b  =  2.5; 

Plot3D[CT,  {tenqjcloudovrbkg,  2,  6},  {f^  -2,  -5}, 

PlotL^l  ->  "Flare  Contrast  Transmittance  Tt~Tb,OD=2.5", 

AxesLabel  ->  {"TfTb",  "log  rfR",  "CT  PlotPoints  ->  25] 
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Program  Darpaflare2 

(*ir  flare  emission  within  a  cloud  is  represented  with  a  high  ten5)erature  \ 
high  albedo  cloud  where  1  - 
albedo  is  the  absorption/ 

emission  which  is  just  the  volume  fraction  of  the  flare  /  (flare  + 
white  smoke)  material  multiplied  by  the  ratio  of  flare  absorption  \ 
coefficient  divided  by  smoke  extinction  coefficient  ~ 
radius  smoke  particles  /  Radius  flare  particles*) 
g  =  0.75; 
p  =  a(l  +  g)/2; 
q  =  a(l  -  g)/2; 
z  =  Sqrt[(l  -  a)(l  -  a  g)]; 

Tu  =  Exp[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Exp[z  b]  -  Exp[-z  b])/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

(♦Note  that  contrast  transmittance  is  just  contrast  ratio  - 
contrast  with  cloud/without  cloud*) 

CT  =  Tu(ten:^targovrbkgM  +  l)/(Tu(terr5)targovrbkgM  +  1)  +  2(R  +  T  -  Tu)  + 
2(1  -  R  -  T)ten:q5cloudovrbkg'^4  ); 

SNT  =  (TuCT)^5; 
b=l; 

a=l-10^ff; 

temptargovrbkg  =  1 .000 1 ; 

Plot3D[SNT,  {tenqicloudovrbkg,  2,  6},  {ff  -2,  -5}, 

PlotLabel  ->  "Flare  S/N  Transmittance  Tt~Tb,OI>=l", 

AxesLabel  ->  {"TfTb",  "log  rfrR",  "SNT  PlotPoints  ->  25]; 
b=1.5; 

Plot3D[SNT,  {tempcloudovrbkg,  2,  6},  {f^  -2,  -5}, 

PlotLabel  ->  "Flare  S/N  Transmittance  Tt~Tb,OD=1.5", 

AxesLabel  ->  {"TfTb",  "log  rfrR",  "SNT  PlotPoints  ->  25]; 
b  =  2; 

Plot3D[SNT,  {tenq)cloudovrbkg,  2,  6},  -2,  -5}, 

PlotLabel  ->  "Flare  S/N  Transmittance  Tt~Tb,OD=2", 

AxesLabel  ->  {"TfTb",  "log  rfR",  "SNT  "},  PlotPoints  ->  25]; 
b  =  2.5; 

Plot3D[SNT,  {tempcloudovrbkg,  2,  6},  {fif  -2,  -5}, 

PlotLabel  ->  "Flare  S/N  Transmittance  Tt~Tb,Or)=2.5", 

AxesLabel  ->  {"TfTb",  "log  rfR",  "SNT  "},  PlotPoints  ->  25] 
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Program  DarpaflareS 

(*ir  flare  emission  within  a  cloud  is  represented  with  a  high  tenperature  \ 
high  albedo  cloud  where  1  - 
albedo  is  the  absorption/ 

emission  which  is  just  the  volume  fraction  of  the  flare  /  (flare  + 
white  smoke)  material  multiplied  by  the  ratio  of  flare  absorption  \ 
coefficient  divided  by  smoke  extinction  coefficient  ~ 
radius  smoke  particles  /  Radius  flare  particles*) 
g  =  0.75; 
p  =  a(l  +  g)/2; 
q  =  a(l-g)/2; 
z  =  Sqrt[(l  -a)(l  -ag)]; 

Tu  =  E)q)[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Ejq)[z  b]  -  Exp[-z  b])/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

(*Note  that  contrast  transmittance  is  just  contrast  ratio  - 
contrast  with  cloud/without  cloud*) 

CT  =  Tu(tenq)targovrbkgM  +  l)/(Tu(temptargovrbkg^4  +  1)  +  2(R  +  T  -  Tu)  + 
2(1  -  R  -  T)tenipcloudovrbkgM  ); 
ten:5)cloudovrbkg  =  4; 
a=l-l(rff; 

ten^targovrbkg  =  1.0001; 

Plot3D[CT,  {b,  1, 2.5},  {fif  -2,  -5}, 

PlotLabel  ->  "Flare  Contrast  Transmittance  Tt~Tb,TfTb=4", 

AxesLabel  ->  {"OD",  "log  rfrR",  "CT  "},  PlotPoints  ->  25]; 
ten5)cloudovrbkg  =  5; 

Plot3D[CT,  {b,  1,  2.5},  {ffi  -2,  -5}, 

PlotLabel  ->  "Flare  Contrast  Transmittance  Tt~Tb,TfTb=5", 

AxesLabel  ->  {"OD",  "log  rfrR",  "CT  "},  PlotPoints  ->  25]; 
tenpcloudovrbkg  =  6; 

Plot3D[CT,  {b,  1, 2.5},  {fif  -2,  -5}, 

PlotLabel  ->  "Flare  Contrast  Transmittance  Tt~Tb,TfTb=6", 

AxesLabel  ->  {"OD",  "log  rfrR",  "CT  "},  PlotPoints  ->  25]; 
ten5)cloudovrbkg  =  3; 

Plot3D[CT,  {b,  1, 2.5},  {ffi  -2,  -5}, 

PlotLabel  ->  "Flare  Contrast  Transmittance  Tt~Tb,TfTb=3", 

AxesLabel  ->  {"OD",  "log  rfrR",  "CT  "},  PlotPoints  ->  25] 
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Program  Darpaflare4 

(*ir  flare  emission  within  a  cloud  is  represented  with  a  high  ten:5)erature  \ 
hi^  albedo  cloud  where  1  - 
albedo  is  the  absorption/ 

emission  which  is  just  the  volume  fraction  of  the  flare  /  (flare  + 
white  smoke)  material  multiplied  by  the  ratio  of  flare  absorption  \ 
coefficient  divided  by  smoke  extinction  coefficient  ~ 
radius  smoke  particles  /  Radius  flare  particles*) 
g  =  0.75; 
p  =  a(l  +g)/2; 
q  =  a(l-g)/2; 
z=Sqrt[(l-a)(l  -ag)]; 

Tu  =  Exp[-b]; 

T  =  2z/((l  -  p  +  z)E;q)[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Exp[z  b]  -  Exp[-z  b])/((l  -  p  +  z)Ejqp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

(♦Note  that  contrast  transmittance  is  just  contrast  ratio  - 
contrast  with  cloud/without  cloud*) 

CT  =  Tu(ten:^targovrbkgM  +  l)/(Tu(ten:5)targovrbkgM  +  1)  +  2(R  +  T  -  Tu)  + 
2(1  -  R  -  T)ten:5)cloudovrbkgM ); 

SNT  =  (TuCT)^5; 
ten:5)cloudovrbkg  =  4; 
a=  1  -  lO^ff; 

ten:5)targovrbkg  =  1.0001; 

Plot3D[SNT,  {b,  1, 2.5},  {fE,  -2,  -5}, 

PlotLabel  ->  "Flare  S/N  Transmittance  Tt~Tb,TfTb=4", 

AxesLabel  ->  {"OD",  "log  rfrR",  "SNT  "},  PlotPoints  ->  25]; 
tencpcloudovrbkg  =  5; 
a=l-10^fiF; 

teraptargovrbkg  =  1.0001; 

Plot3D[SNT,  {b,  1,  2.5),  {S,  -2,  -5}, 

PlotLabel  ->  "Flare  S/N  Transmittance  Tt~Tb,TfTb=5", 

AxesLabel  ->  {"OD",  "log  rfR",  "SNT  "},  PlotPoints  ->  25]; 
ten:5)cloudovrbkg  =  6; 
a=l  -  KTfif; 

ten5)targovrbkg  =  1.0001; 

Plot3D[SNT,  {b,  1, 2.5},  {fE,  -2,  -5}, 

PlotLabel  ->  "Flare  S/N  Transmittance  Tt~Tb,TfTb==6", 

AxesLabel  ->  {"OD",  "log  rfrR",  "SNT  "},  PlotPoints  ->  25]; 
tenpcloudovrbkg  =  3; 
a  =  1  -  lO^ff; 

tenqjtargovrbkg  =  1.0001; 

Plot3D[SNT,  {b,  1, 2.5},  {f^  -2,  -5}, 

PlotLabel  ->  "Flare  S/N  Transmittance  Tt~Tb,TfTb=3", 

AxesLabel  ->  {"OD",  "logrfR",  "SNT  "},  PlotPoints  ->  25] 
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Program  DarpaflareS 

(*ir  flare  emission  within  a  cloud  is  represented  with  a  high  ten:5)erature  \ 
high  albedo  cloud  where  1  - 
albedo  is  the  absorption/ 

emission  which  is  just  the  volume  fraction  of  the  flare  /  (flare  + 
white  smoke)  material  multiplied  by  the  ratio  of  flare  absorption  \ 
coefficient  divided  by  smoke  extinction  coefficient  ~ 
radius  smoke  particles  /  Radius  flare  particles*) 
g  =  0.75; 
p  =  a(l  +  g)/2; 
q  =  a(l  -  g)/2; 
z  =  Sqrt[(l-a)(l-ag)]; 

Tu  =  Exp[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

R  =  q(Exp[z  b]  -  E)q)[-z  b])/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

(*Note  that  contrast  transmittance  is  just  contrast  ratio  - 
contrast  with  cloud/without  cloud*) 

CT  =  Tu(ten:5)targovrbkgM  +  l)/(Tu(ten:5)targovrbkg^4  +  1)  +  2(R  +  T  -  Tu)  + 
2(1  -  R  -  T)ten:5)cloudovrbkg^4  ); 

ff=-2; 

a=l-10^ff; 

tenoptargovrbkg  =  1.0001; 

Plot3D[CT,  {ten:5)cloudovrbkg,  2,  6},  {b,  1,  2.5}, 

PlotLabel  ->  "Flare  Contrast  Transmittance  Tt~Tb,rfrR=.01", 

AxesLabel  ->  {"TfTb",  "OD",  "CT  "},  PlotPoints  ->  25]; 
ff=-3; 
a=l-10^ff; 

ten^jtargovrbkg  =  1.0001; 

Plot3D[CT,  {ten^cloudovrbkg,  2,  6},  {b,  1, 2.5}, 

PlotLabel  ->  "Flare  Contrast  Transmittance  Tt~Tb,rfrR=.001", 

AxesLabel  ->  {"TrTb",  "OD",  "CT  "},  PlotPoints  ->  25] 
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Program  Darpaflare6 

(*ir  flare  emission  within  a  cloud  is  represented  with  a  high  ten:5)erature  \ 
high  albedo  cloud  where  1  - 
albedo  is  the  absorption/ 

emission  which  is  just  the  volume  fraction  of  the  flare  /  (flare  + 
white  smoke)  material  multiplied  by  the  ratio  of  flare  absorption  \ 
coefficient  divided  by  smoke  extinction  coefficient  ~ 
radius  smoke  particles  /  Radius  flare  particles*) 
g=0.75; 

p  =  a(l  +  g)/2; 
q  =  a(l  -  g)/2; 
z=Sqrt[(l-a)(l-ag)]; 

Tu  =  Ejq)[-b]; 

T  =  2z/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Ejq)[-z  b]); 

R  =  q(Exp[z  b]  -  Exp[-z  b])/((l  -  p  +  z)Exp[z  b]  -  (1  -  p  -  z)Exp[-z  b]); 

(♦Note  that  contrast  transmittance  is  just  contrast  ratio  - 
contrast  with  cloud/without  cloud*) 

CT  =  Tu(ten:5)targovrbkgM  +  l)/(Tu(teraptargovrbkgM  +  1)  +  2(R  +  T  -  Tu)  + 
2(1  -  R  -  T)tenq3cloudovrbkgM  ); 

SNT  =  (TuCT)^5; 

ff=-2; 

a=l  -  lO^ff; 

ten:5)targovrbkg  =  1.0001; 

Plot3D[SNT,  {ten:5)cloudovrbkg,  2,  6},  {b,  1, 2.5}, 

PlotLabel  ->  "Flare  S/N  Transmittance  Tt~Tb,rfR.=.01", 

AxesLabel  ->  {"TfTb",  "OD",  "SNT  "},  PlotPoints  ->  25]; 
ff=-3; 
a=l-10'"ff; 

ten^^targovrbkg  =  1.0001; 

Plot3D[SNT,  {tempcloudovrbkg,  2,  6},  {b,  1, 2.5}, 

PlotLabel  ->  "Flare  S/N  Transmittance  Tt~Tb,rfR^=.001", 

AxesLabel  ->  {"TfTb",  "OD",  "SNT  "},  PlotPoints  ->  25] 
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